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Abstract 

We investigate the role of the delta isobar in the reaction ird —> ird at threshold in chiral 
effective field theory. We discuss the corresponding power counting and argue that this calculation 
completes the evaluation of diagrams up to the order x 3 ^ 2 ? with x the ratio of the pion to the 
nucleon mass. The net effect of all delta contributions at this order to the pion-deuteron scattering 
length is 6a£ d = (2.4 ± 0.4) x 1CT 3 m" 1 . 

1 Introduction 

Chiral perturbation theory (ChPT) is the effective field theory of the Standard Model at low energies 
allowing for high accuracy calculations of hadronic observables. It is a systematic expansion around 
the chiral limit (vanishing quark masses) and vanishing external momenta. ChPT can be applied to 
systems containing pions, nucleons and external sources. Here we focus on the ttNN system — in 
particular the pion-deuteron system at threshold. 

In the original formulation, only pions and nucleons appear as dynamical degrees of freedom [H[2], 
whereas the impact of baryon resonances as well as heavier mesons is absorbed into certain low-energy 
constants. From phenomenological studies it is well known that the delta isobar A(1232) plays a very 
special role in low energy nuclear dynamics [3] as a consequence of the relatively large 7riVA coupling 
and the quite small delta-nucleon mass difference A = Ma — ~ 3/^, where Ma, Mjv, and f n 
denote the mass of the delta, of the nucleon, and the pion decay constant, respectively!^! In the effective 
field theory sketched above this leads to unnaturally large values of some low-energy constants. 

It is also possible to include the delta as dynamical degree of freedom in the effective field theory [U 
[5]. For the ttN system this leads to a somewhat improved convergence of the chiral expansion [6], 
however, no qualitative difference appears compared to ChPT. In many cases, the representation of 
delta effects through local pion-nucleon operators is quite accurate. As an example we mention the 

'Tike the pion decay constant, the delta-nucleon mass splitting does not vanish in the chiral limit and thus this 
identification is more appropriate than A ~ 2m„ with m w the pion mass, as often found in the literature. 
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a) b) c) 

Figure 1: Classes of one-body diagrams that contribute to ird scattering. Diagrams a) and b) represent 
the one-body operators with the delta, diagram c) shows the corresponding contact interaction in the 
delta- less theory. Dashed lines denote pions and single (double) solid lines denote nucleons (deltas). 
Solid black dots stand for interactions, whereas the hatched area shows the deuteron wave function. 



successful analysis of threshold pion photoproduction [7j. However, not only for energies of the order 
of A, but also at low energies in the spin sector the explicit inclusion of the delta appears mandatory 

- see, e.g., the recent review [8]. It is, however, important to stress that the delta-full theory in the 
single nucleon sector features more counter terms at a given order than ChPT and has been much less 
systematically applied to low-energy reactions. 

In the present paper we investigate the role of the delta isobar in the reaction ird — > ird at threshold 
in chiral effective field theory. The reason why the explicit inclusion of the delta in pion reactions on 
the two-nucleon system is beneficial is that the irN amplitudes appear in the boosted frame due to 
the Fermi motion of the nucleons. Let us for simplicity focus on the one-body terms with the delta 
that contribute to ird scattering at threshold (see Fig. [1^) and b)). Then the irNA transition vertex is 
linear in the nucleon momentum, p, and the corresponding embedded irN — > irN transition potential 
is proportional to p 2 times the nucleon-delta propagator. The latter behaves as l/(m T — A — p 2 / Mjy) 

— we point out that the width of the delta is suppressed by two powers in the pion mass and thus does 
not contribute to the order we are working. For static deltas, this propagator reduces to l/(m n — A) 
and the sum of diagrams a) and b) of Fig. [1] collapses to diagram c). Thus, in the latter case the 
transition operator behaves like p 2 , whereas in the former it approaches a constant for momenta larger 
than \pj\ I ~ \/(A — m„-)Mjv ~ l.lm^ with the effect that the static amplitude is more sensitive to 
the short range part of the deuteron wave function and must be balanced by appropriate counter 
terms, eventually of unnatural size. The operator with the dynamical delta, on the other hand, does 
not share this problematic property. This point will be discussed in detail below. The value of £>a is 
numerically very close to pthr = \/M^m^ ~ 2.6m,,- — the minimum initial momentum for the reaction 
NN — » NNir. Therefore, in what follows we will use 

PA ~ Pthr > m n . (1) 

It was shown in Ref. [9] that the so-called dispersive corrections to the ird scattering length are 
suppressed by a factor x 3 ^ 2 relative to the leading two-body operator with two Weinberg-Tomozawa 
(WT) vertices, where x = m^/M^. The corresponding power counting, confirmed numerically, treated 
explicitly the scale pthr 3> m n in line with the counting rules for NN — > NNir |10j . The counting rule 
Eq. dU) automatically puts the delta contributions in the same order as the dispersive corrections, as 
we demonstrate below. There is one more class of contributions that can scale as y/x in few-nucleon 
systems, namely the effect of irNN cuts. However, their impact on ird scattering is negligible as shown 
in Ref. [11] . Thus, with this paper we complete the calculation of diagrams at order x 3 / 2 . 
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a) b) c) 

Figure 2: Typical subamplitudes that contribute to irNN scattering including deltas at order x 3 ^ 2 - 



The paper is structured as follows: in the next section we describe the power counting. Results for 
the ird scattering length are reported and compared to previous works in section [3l The paper ends 
with some concluding remarks. 



2 Power counting 

First of all we would like to remind the reader that the leading order two-body operator with two 
WT vertices scales simply as m^/(f^p 2 ) ~ V/Jr in the Weinberg counting scheme, where for this 
diagram p ~ m n . Below we will follow the logic of Ref. [9] and compare diagrams with the delta with 
this leading amplitude. Let us start with the one-body terms depicted in Fig. [Th). The transition 
amplitude scales as 
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As outlined in the Introduction, due to the squared momentum in the numerator, momenta of order 

Ipa 



m. 



-)Mtv contribute to the full matrix element. Therefore the nucleon before (after) 



the pion absorption (emission) is off its mass shell by —p^/2M^ , i.e. by about m w . On the other 
hand, only on-shell amplitudes are physically meaningful and should be compared to each other. To 
find the corresponding chiral order we should therefore estimate the one loop diagram as shown in 
Fig. E^jl. The estimate for this diagram gives 
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where the factors stand for the quantitative estimates for the two one-pion exchange potentials, the 
two two-nucleon propagators, the 7riV — > ttN transition potential through the delta, as defined in 
Eq. ([2]), and the integral measure, consecutively. We stress that we do the power counting based on 
the expressions for time-ordered perturbation theory, since we later work within this formalism. For 
details on this we refer to Appendix E of Ref. [12] . In the relation ([3]) we estimated the contribution 
of diagram [2h.) for two regimes of pion momenta, namely p ~ and p ~ pa- For the identification 
of the chiral order we used Mm ~ 47r/ w . We thus conclude that the power counting yields that the 
dominant contribution of the delta loops is expected to come from loop momenta of the order of p&, as 



2 Note that the external nucleons in Fig. [2^,) can also be off-shell, when the transition operators are convoluted with 
the external wave functions. However, it is the central assumption of the power counting that the corresponding matrix 
element is dominated by (near) on-shell kinematics for these nucleons. 
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a) 




c) 



Figure 3: Typical subamplitudes that contribute to nNN scattering including deltas at order x 2 ■ 
These (and all the others of such kind) are not included in this work. 



argued above. Therefore the considered delta diagram contributes to the same order as the so-called 
dispersive corrections to the ird scattering length as discussed in Ref. [9]. 

The power counting for the other ird diagrams goes in just the same way. E.g., for the diagrams 
b) and c) of Fig. [2] we also find the chiral order x 3 ^ 2 - The power counting for a diagram of type b), 
however with a two-nucleon intermediate state, was discussed in detail in Ref. [Sj. Since we count 
the nucleon-delta propagator in the same manner as the two-nucleon propagator, i.e. as l/m n , it 
becomes obvious that the corresponding diagrams contribute at the same order. In the estimation 
of the chiral order of diagrams b) and c) we used that the leading ttN — > irN vertex scales as m n 
and not as the pion energy, regardless of the relatively large momentum running in the loop. The 
terms dropped are higher order in the chiral expansion. This is in line with the findings of Ref. [13] 
for the reaction NN — > dir. Every additional loop including deltas leads at least to an additional 
factor p 3 /(47r) 2 x l/m w x 1/f 2 for the integral measure, the NA propagator and the leading NA 
interaction, consecutively. Therefore, diagrams with an intermediate NA — ► NA transition, as shown 
in Fig. [3h.) and b), and those with an intermediate A^A" — ► NA transition, diagram c), are suppressed 
by one power in p/Mjy ~ X 1 ^ 2 compared to the diagrams shown in Fig. [2] and will not be considered 
in this work. Consequently, from the naive power counting arguments we can expect the leading delta 
contribution to be of order of (m^/M^) 3 ^ 2 |a^ ublc | ~ 0.06 |a^ p | — 1.6 x 10 _3 m~ 1 where we used that 
|a°^ p | ~ 26 x 10 _3 m~ 1 and that the real part of the scattering length is dominated by the double 
rescattering term with two WT vertices — giving rise to a^ uble . This estimation is fully in line with 
our numerical results as given in the next section. In addition, as stated already, we do not consider 
terms of order x 2 ■ Using the same reasoning as above, we can also estimate the theoretical uncertainty 
of our calculation as 

Aa theor = (m w /M N ) 2 |a^ ublc | ~ 0.6 x lO^m^ 1 . (4) 

At order x 2 also the leading NNtt — > NNtt counter term contributes to nd scattering with up-to-now 
unknown coefficient. Therefore Aa theor represents at the same time an estimate for the theoretical 
accuracy for the extraction of the isoscalar scattering length a + from ird scattering |14| . For a further 
improvement in the accuracy of the calculation, input from other reactions is needed to fix the value 
of the counter term. One possible source of this information could be the reaction AW — ► NNir°ir°. 

3 Results and comparison to previous works 

Although the vertex structure we use for the 7rAA vertex is standard (see, e.g., Ref. [15] and references 
therein) , in order to keep the paper self-contained and to fix the normalization we present it here (note 
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Table 1: Delta contributions to the real part of a^d in units of vn~ l x calculated with Ha = 2.77. 
The results are shown for the phenomenological NN potentials Paris [21], AV18 [22], CD-Bonn [23J 
(CDB), and CCF [23] as well as for the wave functions of the N 2 LO chiral AW interaction [25] based on 
the pairs of regulators {600,500} (EGM1), {550,600} (EGM2) and {450,700} (EGM3). All integrals 
are evaluated up to 1 GeV - the contributions of higher momenta are negligible. 







Paris 


AV18 


CDB 


CCF 


EGM1 


EGM2 


EGM3 


el 


c y 




+1.89 


+1.92 


+1.77 


+1.81 


+1.69 


+1.66 


+1.57 


e2 






+0.54 


+0.55 


+0.73 


+0.56 


+0.72 


+0.79 


+0.84 






















e3 




= 


-0.70 


-0.73 


-0.94 


-0.72 


-0.91 


-1.02 


-1.08 






sum of this group 


= 


+1.73 


+1.74 


+1.56 


+1.65 


+1.50 


+1.43 


+1.33 


/I 

/2 

/3 






+0.84 
+0.13 

-0.05 


+0.85 
+0.14 

-0.05 


+0.75 
+0.21 

-0.14 


+0.79 
+0.14 

-0.06 


+0.69 
+0.21 

-0.14 


+0.67 
+0.24 

-0.19 


+0.63 
+0.26 

-0.22 


< > 








< > 




sum of this group 




+0.92 


+0.94 


+0.82 


+0.87 


+0.76 


+0.72 


+0.67 




total sum 




+2.65 


+2.68 


+2.38 


+2.52 


+2.26 


+2.15 


+2.00 
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that our vertex normalization differs by a factor of two compared e.g. to the one of Ref. [16J): 



J IT 

= _^_[iATt T .^5.v^ A+ h.c.] . (5) 

Here Ha denotes the leading AiW coupling, and S and T are the spin and isospin transition matrices, 
normalized such that 

SiSj = -(25ij - ie ijk a k ) , 

TiTj = -{28ij - ie ijk T k ) . (6) 

In our calculations we use f n = 92.4 MeV and Ha = 3gA/V% — 2-lgA = 2.77, where gA = 1-32 
is the axial-vector coupling of the nucleon (derived from the Goldberger-Treiman relation). The 
relation between Ha and qa can be derived from large N c arguments and the resulting coupling gives 
a reasonable description of the delta width at tree level [17]. Very similar values were shown to be 
consistent with the irN phase shifts in the delta region [18] 119]. It should be noted, however, that the 
dispersion theoretical analysis of Ref. [20J leads to the considerably lower value of Ha = 2.1. 

In Table Q] we show the results of our numerical calculations for the complete set of diagrams with 
the delta isobar that contribute at order x 3 ^ 2 - These numbers were produced using our prefered value 
Ha = 2.77. In order to study the model dependence of the results we performed the calculations for 
various NN potentials. Note that we used phenomenological NN models without [211122^ 123] and with 
[24J explicit delta degree of freedom, as well as three variants of NN wave functions derived within 
chiral effective field theory [25] . We remark that ideally one would also use chiral wave functions with 
explicit deltas. However, up to now corresponding wave functions of sufficient accuracy exist only for 
higher partial waves [16J. Using the different potentials mentioned, we obtain 

5a£ d = (2.38 ± 0.40) x 10~ 3 m" 1 , (7) 

where the central value is the arithmetic average of the results for the seven different potentials and 
the uncertainty reflects the variations in the results. Consistency of the power counting demands that 
the dependence on the NN potential used does not exceed the contribution estimated for the leading 
counter term, Aa theor given in Eq. ([4]), that can absorb this dependence. In this sense Eq. (J7J) is an 
additional confirmation for the consistency of the power counting employed. 

All diagrams evaluated contain the ttNA coupling constant Ha squared. Thus, to see the impact 
of a value as low as Ha = 2.1 on our results, the numbers given in Table Q] simply need to be rescaled. 
We then would get ^^d\h =2 1 = (1-4 ± 0.2) x 10~ 3 m^ 1 . However, we regard Eq. © as our main 
result, since the value employed for Ha can be extracted from fits to the ttN system in the delta region 
based on calculations consistent with the one discussed here [181 ttH] • 

Note that the results from the chiral wave functions are systematically lower than those from the 
phenomenological potentials, which might be a consequence of differences of the NN interactions at 
intermediate range. This finding does not come unexpected. However, calculations to higher orders 
are necessary to draw more firm conclusions. 

The results we found depend only very weakly on the NN models used. In contrast to this, many 
previous works find a significant model dependence when using phenomenological parameterizations 
for some of the diagrams discussed above. For example, in Refs. [26l E7J diagram e3 of Table Q] 
was included by replacing the delta propagator and vertices by the phenomenological irN p-wave 
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Table 2: Results for diagrams el and e3 evaluated with the static delta propagator in units of 10 3 m. 
Integrals are evaluated up to 1 GeV. 







Paris 


AV18 


CDB 


CCF 


EGM1 


EGM2 


EGM3 


el 




+3.0 


+3.1 


+2.5 


+2.8 


+2.2 


+2.1 


+1.9 


e3 




-0.3 


-0.4 


-0.8 


-0.4 


-1.0 


-1.3 


-1.4 



amplitude expressed in terms of the p-wave volumes Co and c\ and evaluated in the boosted frame 
(this is called SP-interference term in Ref. |26j). The evaluated matrix element shows a significant 
model dependence for it scales with the deuteron wave function at the origin (for a more detailed 
study of the model dependence of this quantity see Ref. [28]). To illustrate how large the model 
dependence of the corresponding amplitude could be, we give in Table [2] the results for the diagram e3 
calculated with the static A propagator. The results vary by more than a factor of four when different 
NN models are employed. The corresponding results for the diagram el (see Refs. [261 [271 129] for 
the corresponding phenomenological calculations), also given in Table differ by a factor 1.6. As 
stressed already in the Introduction, once the kinetic energy of the delta is kept in the propagator, 
as demanded by the power counting, the above problem disappears and almost model-independent 
results emerge (see lines el and e3 in Table [I]). In an effective field theory calculation without explicit 
deltas, diagrams el and /l were included effectively as so-called boost corrections [30]. The resulting 
contribution to the 7rd scattering length turned out to be quite sizable, namely (3 — 5) x 10~ 3 m~ , 
depending on the regulator used for the AW potential. Evidently, the spread in the results is well 
above the estimate of Eq. (J3]), which, again, is a consequence of dropping the kinetic energy of the 
delta isobar. In the theory without deltas the pertinent one-body operator scales with the square 
of the nucleon momentum and therefore the corresponding expectation value is proportional to the 
nucleon kinetic energy inside the deuteron — this quantity is strongly model-dependent [28], which 
indicates that the power counting in the delta-less theory requires further modification. However, 
the boost term (see Ref. [30]) is proportional to the low energy constant C2, which is known to be 
largely saturated by the delta isobar [31[ 116] . In the analysis of the ttN system [31] it was shown 
that the explicit evaluation of the leading order delta contribution results in a reduction of the value 
of C2 from about 3.3 GeV -1 to about 0.5 GeV -1 . In the very recent analysis of the NN system 
including explicitly the delta at NLO |16j . an analysis of ttN threshold coefficients was performed. 
Given the parameters utilized there, the value of C2 is reduced to —0.25 GeV -1 . A reduction of the 
irNA coupling by 30% as demanded by a dispersive analysis of the resonance contribution to the 
pion-nucleon P33 phase shifts [20] leads to a reduced C2 = 0.83 GeV -1 . All these values are consistent 
within the uncertainty of the various contributions to the low-energy constants given in Ref. [31] , 
Therefore the value of C2 is reduced by a large factor once the delta contribution is taken out. We 
have calculated the residual boost correction using the expressions given in Ref. [3D] with N 2 LO chiral 
wave functions and with C2 = —0.25 GeV -1 and found it to be as small as —(5.7. . . 6.6) x 10~ 4 m~ . 
Consequently, this correction is of the same size as the estimated uncertainty of the calculation (see. 
Eq. ([!])) and thus does not contribute significantly anymore. 
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It should be stressed that it is not compulsory for a consistent calculation of the nd scattering length 
that the delta is included explicitly. Also a calculation without deltas is obviously equally justified. 
As usual the effects of the delta would then be parameterized by local counter terms of the type 
ttNN — > ttNN with up-to-now unknown coefficients. The conclusion to be drawn from our studies is 
that in order to perform calculations with the accuracy of the order of the uncertainty estimate given 
in Eq. ([4]) it is necessary to include a dynamical delta, as long as no additional information on the 
size of the counter term is available. On the other hand, for a consistent inclusion of isospin breaking 
effects, that are known to be important [32J, more theoretical work on the treatment of effects from 
quark masses and virtual photons in the delta-full theory would be useful. 

4 Conclusions 

In this work we calculated the leading contributions of the A(1232) to the ird scattering length in 
effective field theory. As expected, inclusion of the delta leads to an improved convergence for the 
isospin-symmetric operators that contribute to this reaction. We have also compared our results to 
other approaches and discussed the differences. 

In the power counting employed the delta starts to contribute at order x 3//2 > relative to the leading 
two-nucleon contribution, given by two subsequent ttN scatterings on the two different nucleons. At 
the same chiral order the so-called dispersive corrections evaluated in Ref. [9] contribute as well, 
and with this work we complete the evaluation of diagrams at that order. In Ref. [9] the dispersive 
corrections were evaluated for a particular NN potential. When repeating the calculation with the 
four different phenomenological NN potentials employed in the present study (note: the chiral wave 
functions could not be used here, since for the dispersive corrections the wave functions are needed 
also at pion production threshold, where the chiral wave functions are not applicable anymore) we 
find 

,5^7 = (-2.9 ± 1.4) x lO-Sm" 1 , (8) 

where the first number is the mean value for the various potentials and the second number reflects 
the theoretical uncertainty of this calculation estimated conservatively — see Ref. [9] for details. The 
variation of the results for the different potentials lies well within this uncertainty band. Note, that 
the uncertainty can be reduced by a calculation of NN — > dir to next-to-next-to-leading order, which 
is planned for the near future. We therefore find for the total contribution at order x 3 ^ 2 

6a£ d + 6af° v = (-0.6 ± 1.5) x 10~ 3 m' 1 , (9) 

where we added the uncertainties given in Eqs. (J7|) and (jSJ) in quadrature. Thus, we conclude that 
the net effect of the diagrams that contribute at order x 3 ^ 2 is very small. Note that the occurring 
cancellation is accidental because very different physics contributes to the two classes of diagrams. 

One important consequence of our investigations is that once the delta isobar is treated dynamically, 
as it is done in this paper, the so-called boost corrections give rise to an insignificant contribution in the 
theoretical analysis of the nd scattering length. Furthermore, for the same reason the phenomenological 
inclusion of pion rescattering (the so-called SP interference term) through a boosted p-w&ve amplitude, 
used in Refs. [261 I27j . is expected to yield a very small contribution, well within the theoretical 
uncertainty given here — see also the corresponding discussion in Ref. [9]. 

With this work all strong, isospin-symmetric contributions to the ird scattering length have been 
calculated to very high accuracy. In principle we could now extract the isoscalar ttN scattering length, 
a + , directly from the ird scattering length, since 

a n - d = 2a + + (few-body corrections (a - )) , (10) 
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where a~ denotes the isovector scattering length. In this expression additional terms that contain a + 
were neglected for they are numerically negligible. However, in addition isospin violating effects are 
known to be quite sizable. Therefore, in Eq. ([TO]) we should replace 2a + by a n -„ + a n ~ n which agrees 
to the former only, if isospin were an exact symmetry. Furthermore, few-body corrections involving 
virtual photons, in addition to those calculated in Ref. [9], are potentially important. For the n~d 
system so far only the leading isospin violating corrections were evaluated [32] . To this order the 
largest theoretical uncertainty emerged from the appearance of the low-energy constants f\ and c\. It 
is intriguing to observe, however, that those appear in the same linear combination in both a w - p and 
a^- n . Thus, one is in the position to extract a n - p + a n - n with high accuracy from a combined analysis 
of pionic deuterium and pionic hydrogen even without detailed knowledge on /i (see also Ref. [33]). 
However, it remains to be seen if the corrections at next-to-leading order in isospin violation do not 
distort this picture. Corrections at this order for the n~p system were evaluated in Refs. [34^ [35] and 
turned out to be quite sizable, especially those that come from the pion mass difference. In order 
to push also the calculation for the nd system to a similar level of accuracy in isospin violation, the 
ir~n scattering amplitude as well as some virtual photon exchanges in the ir~d system are still to be 
calculated. 
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Appendix 



In this appendix we present the explicit expressions for the amplitudes given in Table [TJ Note that 
in accordance with the power counting, we only keep those amplitudes that contain intermediate 
states with the nucleon, the delta and at most real pions. The calculation is done in time-ordered 
perturbation theory (TOPT). Especially, we dropped the so-called stretched boxes. The corresponding 
correction to the ird scattering length due to the delta isobar is 



Sa 



(Al) 



where the first and second terms correspond to the direct and crossed diagrams of Tabled! respectively. 
Here 



a nd(Q0) 



h\ml 



d\ 



q- 



(h Pl +h Pl +h Fl ) (A.2) 



48^/6 (l+m w /2M N ) J (2tt)* qo - A - q z /2M NA 

where Mjva = M^M^/ (Mn + Ma) and hs+i Lj are the partial wave amplitudes squared that corre- 
spond to the decomposed intermediate NA state 



I 



3 Pi 



l5 Fl 



2V2 
3 



lOv^ 
2f? 



4(q) - 



M A 



u(q) + 
u(q) - 



w{q) 

V2 

w(q) \ 
' 5V2j 



n 2 



(A.3) 



M A 



w{q) 
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with the 1^ denoting the integrals that correspond to the overlap of the deuteron wave function (u(q) 
and w(q) for the S- and D-waves, respectively) with the one-pion-exchange operator 

T A, , f <?p l-{p-q?/{p 2 q 2 ) . . / 1 1 \ 

A /■ d 3 p 3 + 4(p-g)/g 2 + (p-g)V(pV) / 1 , 1 
/3((?) = "7(2^ 2^ ^U + if 



A 1 f d*p -1 - 3(p • qW + 3(P • gTAgV) + 5(p • gr/bV) / v 
/4(9) = 2^ 

Here Pi and P^ correspond to the TOPT components of the pion propagator 

Pi = qo -^ +tf -(p 2 + q 2 )/2M N , 

P 2 A = -LO f+(! -A-p 2 /2M N -q 2 /2M A (A.5) 



with ujq- = \Jq 2 + m 2 . The diagrams of Table [T] can be easily matched to the individual terms of 
Eqs. (|A.2p and ()A.3j) : the very last terms on the r.h.s. of each amplitude hs+i Lj in Eqs. ()A.3[) . 
proportional to the deuteron wave functions squared, correspond to the diagrams of type 1 (el and 
/l), type 2 contains If' amplitudes squared, whereas the interference terms of type 3 contain the rest. 
For the direct terms, labeled as e in Table [fl one needs to take qo = and for the crossed terms, 
labeled as / in that Table, qo = —m^. Finally, we remark that all integrals are evaluated up to a 
sharp momentum cut-off of 1 GeV. All higher momentum contributions are negligible and anyway are 
to be absorbed in a counter term that is to be included at order x 2 - calculated with different wave 
functions demonstrates nice convergence. 
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